Observation of (892) "K* (892)° in Xcj Decays 
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K* {892)° K* {892)° signals from XcJ {J ~ 0,1,2) decays are observed for the first time using 
a data sample of 14 million ^p{2S) events accumulated at the BESII detector. The branching 
fractions B{xcj K* (892)° K* (892)°) (J = 0, 1,2) are determined to be (1.55±0.35±0.30) x 10~^ 
(1.58±0.32±0.29) x 10"^ and (4.67±0.55±0.85) x 10"^ for the Xco, Xci and Xc2 decays, respectively, 
where the first errors are statistical and the second are systematic. The significances of these signals 
are about 4.2(t, 4.3cr, and 7.5a, respectively. 

PACS numbers: 13.25.Gv, 12.38. Qk, 14.40.Gx 
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I. INTRODUCTION 

Exclusive quarkonium decays constitute an interesting laboratory for investigating perturbative quantum chromo- 
dynamics (QCD). In the case of P-wave charmonium Xcj decays to a pair of pseudoscalars, one finds that the lowest 
Fock state, the color-singlet contribution, alone is not sufficient to accommodate the data. Indeed, it turns out that 
the color-octet contribution from the next higher Fock state contributes at the same level as the color singlet one. 
Its inclusion yields good agreement with experimental data [1,2]. The calculation of the partial width of Xcj ~^ PP, 
taking into account the color octet mechanism j2i , also obtains results in reasonable agreement with measurements jj] . 
Nevertheless a recent measurement of the Xcj AA |^ only agrees marginally with this prediction. 

At present there are no predictions for the majority of the hadronic decay modes. In addition, few two-body decays 
have been measured. A consistent set of predictions for the branching fractions, as well as more precise experimental 
measurements, for a number of the two-body decays may lead to further insight into the nature of these '^Pj cc bound 
states. 

In this paper, we report on the analysis of Tr+Tr^if+A'^ final states from Xcj {J = 0,1,2) decays using 14 
milhon i{^{2S) events accumulated at the upgraded BES detector (BESII). Signals of XcOj Xa and Xc2 decays to 
K*{892fK*{892)° in i/'(25') radiative decays are observed for the first time. 

II. BES DETECTOR 

BES II is a large solid-angle magnetic spectrometer that is described in detail in Ref. . Charged particle momenta 
are determined with a resolution of cTp/p — 1.78%-\/l -|- {p in GeV/c) in a 40-layer cylindrical drift chamber. Particle 
identification is accomplished by specific ionization (dE/dx) measurements in the drift chamber and time-of- flight 
(TOE) measurements in a barrel-like array of 48 scintillation counters. The dE/dx resolution is (JdEjdx = 8.0%; the 
TOE resolution is cttof — 180 ps for Bhabha events. Outside of the time-of-flight counters is a 12-radiation- length 
barrel shower counter (BSC) comprised of gas proportional tubes interleaved with lead sheets. The BSC measures the 
energies of photons with a resolution of cte/E ~ 21% /^/E {E in GeV). Outside the solenoidal coil, which provides a 
0.4 T magnetic field over the tracking volume, is an iron flux return that is instrumented with three double layers of 
counters that are used to identify muons. 

In this analysis, a GEANT3 based Monte Carlo simulation package (SIMBES) with detailed consideration of detector 
performance (such as dead electronic channels) is used. The consistency between data and Monte Carlo has been 
checked in many high purity physics channels, and the agreement is quite reasonable. 

III. EVENT SELECTION 

The selection criteria described below are similar to those used in a previous BES analysis 0. 

A. Photon identification 

A neutral cluster is considered to be a photon candidate when the angle between the nearest charged track and 
the cluster is greater than 15°, the first hit is in the beginning six radiation lengths, and the difference between the 
angle of the cluster development direction in the BSC and the photon emission direction is less than 30° . The photon 
candidate with the largest energy deposit in the BSC is treated as the photon radiated from ^{2S) and used in a 
four-constraint kinematic fit to the hypothesis '0(25') ^tt'^tt^ . 

B. Charged particle identification 

Each charged track, reconstructed using the MDC information, is required to be well fit to a three-dimensional 
helix, be in the polar angle region | cos 9m dc\ < 0.80, and have the point of closest approach of the track to the beam 
axis be within 2 cm of the beam axis and within 20 cm from the center of the interaction region along the beam line. 
Eor each track, the TOF and dE/dx measurements are used to calculate x^ values and the corresponding confidence 
levels for the hypotheses that the particle is a pion, kaon or proton (Prober, Prob^, Probp). 
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C. Event selection criteria 

Candidate events are required to satisfy the following selection criteria: 

(1) The number of charged tracks is required to be four with net charge zero. 

(2) The sum of the momenta of the two lowest momentum tracks is required to be greater than 650 MeV; this 

removes contamination from ip{2S) tt+tt" J/i/i events and some of the p^tttt background. 

(3) The probability for the four-constraint kinematic fit to the decay hypothesis ip{2S) jTr~^7r~K~^K~ is 
required to be greater than 0.01. 

A combined probability determined from the four-constraint kinematic fit and particle identification information is 



used to separate ■^tt^tt 



',^K^K K^K , and the different possible particle assignments for the 77r"*"7r K^K 



final states. This combined probability, Probow, is defined as 



ProbaH = Prob(Xa/;,ndfa;/), 

where x^;; is the sum of the values from the four-constraint kinematic fit and those from each of the four particle 
identification assignments, and ndfaj; is the corresponding total number of degrees of freedom. For an event to be 
selected, I'mhaii of the "f'K^'K~ K~ must be larger than those of the other possibilities. In addition, the particle 
identification probability of each charged track Prob/75 must be > 0.01. 

The invariant mass distribution for the K^K~ events that survive all the selection requirements is shown in 

Fig. 1. There are clear peaks corresponding to the XcJ states. The highest mass peak corresponds to charged tracks 
final states that are kinematically fit with an unassociated, low energy photon. 




M(7i"'7i"K"'K") (GeV/c^) 
FIG. 1: The ■k^v~ invariant mass spectrum. 



IV. ANALYSIS RESULTS 



The scatter plots of K 7r+ versus K'^-k invariant masses for events with a tt+tt K^K mass within (3.30, 3.48) 
GeV, (3.48, 3.53) GeV and (3.53, 3.65) GeV are shown in Fig. 2. Clear K* K* {S%2f signals can be seen in 
all XcJ decays, as well as some hints of J!'J_(1430)ii'^(1430) (or ii'^(1430)/i'^(1430)) and Ki{l27Q)K (or Kx{lAm)K) 
signals. In this paper, we study ii'*(892)°ii'*(892)° production in Xc0,i,2 decays. 



A. if*(892)°if*(892)° signal 

For the events in Xcj mass region (3.30, 3.65) GeV, after requiring that the mass of either (or both) Kt: pair lies 
between 0.836 and 0.956 GeV, the mass distribution of the other Kir pair, shown in Fig. 3, is obtained; there is a 
strong K* (892) signal. The distribution is fitted with a background polynomial plus a P-wave relativistic Breit-Wigner 
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FIG. 2: Scatter plots of K versus K^-n invariant masses for selected 77r+7r K events with tt+tt K mass in 
XcO, Xci, and Xc2 mass regions, respectively. 



function, with a width 



r = ro 



mo 1 



TO 1 + r2p2 



where m is the mass of the Kt: system, p is the momentmii of kaon in the Kt: system, Fq is the width of the resonance. 

Too is the mass of the resonance, po is P evaluated at the resonance mass, r is the interaction radius, and 

1 + r^p^ 

represents the contribution of the barrier factor. The fit of Fig. 3 gives an r value of (3.4 ± 2.6) GeV~^ with a large 
error due to the low statistics. Therefore, in later analysis (mainly in the efficiency calculation), we use the value 
(3.4 ± 0.6 ± 0.3) GeV^^ measured by the K~it'^ scattering experiment ^ for r. 
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FIG. 3: K-n invariant mass spectrum recoiling against if* (892) (0.836 GeV < nrtK-w < 0.956 GeV) for events in the Xcj mass 
region (some events appear twice), where the curves are the Breit-Wigner function and background polynomial described in 
the text. 



In this paper, the number of if * (892)" A"* (892)° events and the corresponding background are estimated from the 
scatter plot of if ~7r+ versus if +7r^ invariant masses, as shown in Fig. 4. The signal region is shown as a square box 
(solid line) at (0.896, 0.896) GeV with the width of 60 MeV. From a Monte Carlo study, a large background comes 
from il){2S) -fXcj ^ 7ifi(1270)if (or ifi(1400)if) which decays to 77r+7r-/f +/f- final states via ifi if*(892)7r 
intermediate decay. This background shows up as the horizontal and vertical bands at TO(if*(892)) in the TO(if~7r^) 
versus TO(if+7r~) scatter plots of Fig. 2. Hence, backgrounds are estimated from sideband boxes, which are taken 
60 MeV away from the signal box and shown as four dashed-line boxes in Fig. 4. Background in the horizontal or 
vertical sideband boxes is twice that in the signal region. 
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FIG. 4: Definition of signal and sideband regions. 



Figure 5 shows the mass distribution of the if * (892)" A"* (892)° candidate events and the corresponding background. 
There are 154 and 46 (92/2) events obtained from the signal box and the dashed-hne boxes within 3.20 GeV - 3.70 GeV, 
respectively. 

> I • ' • ' • ' • ' • I 
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M(K*(892)°K*(892)") (GeV/c^) 

FIG. 5: The Jf* (892)" A"* (892)" invariant mass spectrum for events in the signal region of Fig. 4. The shaded histogram 
indicates the distribution of the background estimated for events in the sideband regions of Fig. 4. 



B. Fit of the mass spectrum 

After sideband subtraction, the K* (892)'^ K* (892)'^ mass spectrum between 3.20 and 3.70 GeV is fitted using a 
method with three Breit-Wigner functions folded with Gaussian resolutions, where the mass resolutions are fixed at 
their Monte Carlo predicted values [(12.2 ± 0.4) MeV, (12.3 ± 0.3) MeV and (12.2 ± 0.3)MeV for Xco, Xa and Xc2, 
respectively] and the widths of the three Xcj states are set at their world average values 01 • A probability of 70% 
is obtained, indicating a reliable fit. The number of events determined from the fit are 26.1 ± 5.8, 26.9 ± 5.4 and 
55.1 ± 6.3 for XcO, Xci, and Xc2, respectively. The statistical significances of the three states are 4.2fT, 4.3cr and 7.5(J, 
calculated from Ax^, where Ax^ is the difference between the x^ values of the fits determined with and without 
the signal function. Fig. 6 shows the fit result, and the fitted masses are 3416.2 ± 3.6 MeV, 3507.8 ± 3.6 MeV and 
3553.6 ±1.8 MeV for XcO, Xci and Xc2, respectively, in agreement with the world average values 
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M(K*(892) V(892)°) (GeV/c^) 

FIG. 6: The K* (892)" K* {892f invariant mass spectrum fitted with three resolution smeared Breit-Wigner functions as de- 
scribed in the text. 



A Monte Carlo simulation is used to determine the detection efficiency. The angular distribution of the photon 
emitted in ip{2S) iXca is taken into account The if* (892) is generated as a P-wave relativistic Breit-Wigner 
with r as 3.4 GeV^^ 8]. For each case, 50,000 Monte Carlo events are simulated, and the efficiencies are estimated to 

where the error is the statistical error of 



be e^.o = (3.15 ±0.09)%, e^^, = (3.25 ± 0.09)%, and e^.^ = (2.96 ± 0. 
the Monte Carlo sample. Note that for the efficiency estimation, the events in the four sideband boxes are subtracted 
from the events in the signal region of the scatter plot, similar to the treatment of data. 
The branching fraction of 'ip{2S) -fXcj, Xcj ^ if * (892)°^* (892)" is calculated using 



BiXcJ ^ 7XcJ ^ 7i^*(892)°if*(892)") 



4,(23) 



where the factor / is the branching fraction of A'* (892)° to the charged Ktt mode, which is taken as ^. 

Using the numbers obtained above and the total number of tp{2S) events, 14.0 (1.00±0.04) x 10^ '^^ determine 
the following branching fractions 

B{^{2S) 7Xco ^ 7is:*(892)°A:*(892)°) = (1.33 ± 0.30) x 10-^ 

Bi^pi2S) 7Xci ^ 7is:*(892)''A:*(892)°) = (1.33 ± 0.27) x 10"^ 

B{^p{2S) 7Xc2 7is:*(892)''if *(892)°) = (2.99 ± 0.35) x 10-^ 
where the errors are statistical only. 

C. Systematic errors 

The systematic errors in the branching fraction measurement associated with the efficiency are determined by 
comparing ip{2S) data and Monte Carlo simulation for very clean decay channels, such as ip{2S) vr+vr^ J/-0, 
which allows the determination of systematic errors associated with the MDC tracking, kinematic fitting, particle 
identification, and efficiency of the photon ID pdl |. Other sources of systematic error come from the uncertainties 
in the number of ipi'^S) events the efficiency estimation using simulated data, the background, the XcJ and 
if* (892)" mass resolutions, the binning and fit range, etc. 



1. Ejficiency estimation 

As mentioned above, we use the measurement of Ref. 0, (3.4 ± 0.6 ± 0.3) GeV~^ for r in the P-wave relativistic 
Breit-Wigner parameterization in the Monte Carlo simulation. We also use r varied by one sigma to 2.73 GeV~^ and 
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4.07 GeV ^ to determine the change in the detection efficiency. For r = 2.73 GeV ^, the efficiencies of the Xco, Xci, 
and Xc2 become 2.90%, 2.99% and 2.74%, and for r = 4.07 GeV^\ the corresponding efficiencies are 3.37%, 3.52%, 
and 3.16%, respectively. The largest changes are about 7.9%, 8.3% and 7.4% for XcO, Xa and Xc2- 

2. Background subtraction 

In Section IV A, the backgrounds are estimated using the sidebands shown as the four dashed-linc boxes in Fig. 4. 
Moving the sideband boxes 20 MeV away from or closer to the signal region, or varying the background number by 
one standard deviation, the largest changes of the branching fractions for the XcO, Xci and Xc2 are about 7.4%, 5.0% 
and 5.5%, respectively, obtained by re-fitting the (892)" A"* (892)° mass spectrum and reestimating the efficiency. 

3. Xc.J anrf if* (892)" mass resolutions 

Differences between data and Monte Carlo for the mass resolutions of the XcJ or K*{892)'~' also give uncertainties 
in the determination of the branching fractions. The maximum possible difference for Xcj is about 1 MeV. Such a 
change results in about 4.5%, 2.5% and 2.0% variations in the fitted number of XcO: Xci: and Xc2 events. If we change 
the (892)0 window to [0.836 + 0.002, 0.956 - 0.002] GeV and [0.836 - 0.002, 0.956 + 0.002] GeV, the efficiency 
variations of the XcO: Xci and Xc2 are 1.5%, 2.5% and 2.4%, respectively. By varying the width of XcO by la, 0.8 
MeV, there is almost no change in the final fit result. We use total systematic errors of 5%, 3.5% and 3.5% for this 
uncertainty. 

4- Binning and fit range 

Using different binning and fit ranges for the if * (892)" A"* (892)" mass spectrum fit yields errors of about 4%, 2% 
and 2% for XcOi Xci and Xc2, respectively. 

TABLE I: Summary of systematic errors in the branching fraction calculation of B{'tp{2S) ^ 7Xc,7 7^"* (892)"^" (892)"). 



Source 


XcO 


Xci 


Xc2 


MDC tracking 


8% 


8% 


8% 


Kinematic fit 


6% 


6% 


6% 


Particle identification 


5% 


5% 


5% 


Photon ID efiiciency 


2% 


2% 


2% 


ik{2S) number 


4% 


4% 


4% 


Efficiency estimation 


7.9% 


8.3% 


7.4% 


Background 


7.3% 


5.0% 


5.5% 


Mass resolutions 


5% 


3.5% 


3.5% 


Binning and fit range 


4% 


2% 


2% 


Total systematic error 


17.4% 


16.0% 


15.7% 



The systematic errors from all sources are listed in Table I, as are the total errors of 17.4%, 16.0% and 15.7% for 
XcOj Xci and Xc2j respectively, obtained by adding them in quadrature. The resulting branching fractions are 

B{'il^{2S) 7Xco 7is:*(892)°/r(892)") = (1.33 ± 0.30 ± 0.23) x 10~^ 

B{i}^{2S) 7Xci ^ 7if *(892)°if *(892)") = (1.33 ± 0.27 ± 0.21) x 10"^ 

^(^^(25) ^ 7Xc2 ^ 7if *(892)°jr(892)") = (2.99 ± 0.35 ± 0.47) x 10-^ 
and with the PDG world average values of ip{2S) JXcj 4], we get the following branching fractions 

Bixco ^ if *(892)°if *(892)") = (1.55 ± 0.35 ± 0.30) x 10~^ 
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TABLE II: Summary of numbers used in the branching fraction calculation and branching fraction results. 



Quantity 
e(%) 

f 

B(V(2S) jxcj) (%) f4] 

^if*(892)0i?*(892)0) (lO'^) 



XcO 

26.1 ± 5.8 
3.15 ±0.09 
14.0 ±0.6 
2 

8.6 ±0.7 
1.55 ±0.35 ±0.30 



Xcl 

26.9 ± 5.4 
3.25 ±0.09 
14.0 ±0.6 
2 

8.4 ±0.8 
1.58 ±0.32 ±0.29 



X c2 

55.1 ±6.3 
2.96 ± 0.08 
14.0 ± 0.6 
2 

6.4 ±0.6 
4.67 ±0.55 ±0.85 



B{xci K*{892fK*{892)°) = (1.58 ± 0.32 ± 0.29) x 10-^ 

B{xc2 K*{892)°K*{892f) = (4.67 ± 0.55 ± 0.85) x 10-^ 

where the first error is statistical and the second is systematic. The numbers used and results are summarized in 
Table IL 



V. SUMMARY 

In summary, K* (892)'^ K* (892)^ signals from Xd {J = 0,1,2) decays are observed for the first time using a 
sample of 14 million i}{2S) events accumulated at the BESII detector. The branching fractions are determined 
to be B{xcQ -> K*{892.f K*{892f) = (1.55 ± 0.35 (stat) ± 0._30 (syst)) x lO^^, B{xci ^ K*{892f K*{892f) 
= (1.58±0.32 (stat)±0.29 (syst)) x lO^^, and B{xc2 ^ K*{892fK*{892)°) = (4.67±0.55 (stat)±0.85 (syst)) x lO^^. 
and the significances of the K* (892)'^ K* {892)'^ signals are about 4.2(t, 4.3(t and 7.5a for the Xc0,i,2 decays, respectively. 
This will be helpful in understanding the nature of Xc states. 
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